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Abstract. The Kohn-Sham (KS) equations determine, in a self-consistent way, 
the particle density of an interacting fermion system at thermal equilibrium. We 
consider a situation when the KS equations are known to have a unique solution 
at high temperatures and that this solution is a uniform particle density. We 
prove that, at zero temperature, there are stable solutions that are not uniform. 
We provide the general principles behind this phenomenon, namely the conditions 
when it can be observed and how to construct these non-uniform solutions. Two 
concrete examples are provided, including fermions on the sphere which are shown 
to crystallize in a structure that resembles the C60 molecule. 



1. Introduction 

We consider a system of N interacting fermions in a finite volume. Since we want 
to avoid the surface effects, we actually consider the fermions moving on toruses and 
spheres, or, more generally on a closed Riemann manifold M. of finite volume O. 
According to the Kohn-Sham theory (and later extensions 0), the particle density 
at thermal equilibrium at a temperature T (j3 = 1/kT) is a solution of the following 
set of equations: 

{H + V n )4>i = Ci & (1) 
n(z) = ]T(l + e^-^)~W)| 2 , (2) 

i 

with V n an effective potential depending entirely on the particle density n and 

fi determined from J n(x)dx = N. Hq is the single particle, non-interacting 
Hamiltonian. We refer to 

H n = H Q + XV n (3) 

as the Kohn-Sham Hamiltonian, where we introduced the coupling constant A for 
convenience. We will neglect the spin degree of freedom. 

There is no closed form of V n . However, at least for the electron gas, there is a set 
of very successful, explicit approximations, which already provide numerical results 
that are within the so called " chemical precision" |3] . Although in this paper we do not 
use a specific approximation, we will often make reference to two such approximations 
in order to check our assumptions. They are the Local Density Approximation (LDA): 
V n — v * n + v xc (n) (v the two body interaction), i.e. V n is the sum between the 
Hartree potential and a function of density, and the quadratic approximation (QA): 
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V n = K * n, i.e. V n is the convolution of the density with a certain kernel. The 
mathematical structure of QA is the same as that of the Hartree approximation. 

In this paper, we are not concerned with the physical and mathematical principles 
leading to the KS equations, but rather with the mathematical structure of these 
equations, in particular with the question of uniqueness at zero temperature. In 
other words, we have already assumed u-representability, picked our approximation 
for V n and we are ready to compute the self-consistent solutions. What should we 
expect? Well, as already demonstrated (see for example Ref. @]), even for local 
approximations, we should expect a very rich structure, which may include multiple 
solutions, symmetry breaking etc.. While previous studies used purely numerical 
methods, here we use group theoretical methods and functional analysis to study this 
structure. 

Let us discuss first what it is known about the KS equations at finite temperature 
and finite volume. The notation ||</>|| lp stands for [J M \cf>(x)\ p dx] 1 / p and <j) £ LP means 
\\4>\\lp < oo. Assume the following: 

Al) Hq is self adjoint, bounded from below; for —a below its energy spectrum, the 
kernel (H + a)~ 2 (x, x') is continuous (with respect to x and x') and 

k a = sup (H a + a)~ 2 (x, x) < oo. (4) 

A2) V n £ L 2 (Ai) and w = sup ||Kj||j,2 < oo, where the supremum is taken over all n 
in 

S N = {neL 1 (M), \\n\\ L r=N}. (5) 

As Hq is in general equal to minus the Laplace operator, Al is easy to check for 1, 2 
and 3 dimensional toruses or spheres. It fails in 4 and higher dimensions. Since M. 
is of finite volume, Al automatically implies that exp(— (3Hq) is trace class. Al also 
implies that < Vk^\\(Ho + a)f\\ L 2. Together with A2 (easy to verify for LDA 

and QA jSJEj), this leads to 

\\V n (H +a)- 1 \\<wVk~a = la. (6) 

Then, H n is self-adjoint for all n £ S and, as it follows from Ref. 0, the Kohn-Sham 
equations can be formulated as a fixed point problem: 

Theorem 1. For T > 0, the following map is well defined: 

T : S N — > S N 

S N 3 n ^ T[n]{x) = (l + e ^ H "-^ ^ (x, x), (7) 

where fj, is the unique solution of N — Tr (l + e /3< - Hn ~^) 1 . The fixed points of T 
generate all possible solutions of the KS equations. 

Many will recognize in Eq. (0 the usual formulation of the KS problem in terms 
of the density matrix. When appealing to the fixed point theorem, the functional form 
of the map T and its domain of definition are equally important. What is new in the 
above result is that T is well defined for all densities which integrate to N. 

Apart from complications that may occur at low particle densities and which will 
not be addressed here, the following assumption can be easily verified for LDA and 
QA (see Ref. 0EJ): 
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A3) There exists x < 00 such that 

\\V n -V n ,\\ L t< X \\n-n'\\ L t, (8) 
for any n, n 1 E S . 

If A1-A3 are satisfied, then there exists k, which is a function of A, such that [5] 

\\T[n]-T[n']\\ Ll <K\\n-n'\\ L r, (9) 

on the entire S . For A smaller than a critical value A c , T becomes a contraction and, 
consequently, it has a unique fixed point. If the constants can be chosen independent 
of temperature in A1-A3, it is not hard to show that A c increases with temperature. 
In other words, if A is kept fixed, A1-A3 (and the fact that f2 is finite) guaranties the 
existence of a unique fixed point of T at high temperatures. 

Let us end the finite temperature case with a few remarks. For an exact V n , 
the existence and uniqueness, at any finite T, will follow from the convexity of the 
functional [7], provided that the equilibrium density can be written as in Eq. @ (i.e. 
is w-representable) . In practice, we don't have the exact V n and v-representability has 
not been yet proved or disproved. Also, in the thermodynamic limit, where systems 
can have multiple coexisting phases, the issue of uniqueness becomes more delicate 
and definitely there are many opened questions here. Thus, the question of existence 
and uniqueness in the finite temperature Kohn-Sham equations is not trivial. 

The situation at zero temperature is more delicate. The density now becomes 
n(x) = ^2 \4>i( x )\ 2 where the sum goes over the lowest N energy states of H n . If 
the last occupied energy level is degenerate and only partially occupied, there is an 
ambiguity in defining n(x). In this paper, we deal exactly with this situation. 

Let us assume that that there is a continuous group G c acting ergodically on M. 
and preserving the Riemann structure. On torus or sphere, this group will be simply 
the translations or rotations. Let us consider the natural unitary representation of G c 
in L 2 (M): 

G c 3 9^g, (gf)(x) = f(gx). (10) 

We assume that Hq commutes with all g and that every symmetry of the particle 
density is automatically a symmetry of the effective potential: 

A4) If n(gx) = n(x), then V n (gx) = V n (x) (equivalently gV n g^ = V n ). 

This assumption can be easily verified for LDA and QA. Besides other things, A4 
implies that V n is a constant if n{x) is uniform, and we can fix this constant to zero. 
In other words, the Kohn-Sham Hamiltonian reduces to Hq if n(x) — h (h — N/tt). 
Then, it is trivial to show that, at any finite temperature, n is a solution of the KS 
equations. With our assumptions, we also know that this is the only solution at high 
temperatures. At zero temperature, assume that, if we populate with N particles the 
energy levels of Hq, from smaller to higher energies, we end up with Nq particles on the 
last occupied energy level, assumed d-fold degenerate with d > N . We refer to this 
level and its energy as the Fermi level and Fermi energy e® F . If we can find Nq states 
at the Fermi level so as to generate a uniform particle density, then n is a solution 
of the KS equations. If there is no such combination of states, than either there is 
no solution, the solution is not uniform or we need to consider fractional occupation 
numbers. We will not discuss here the last possibility, but rather concentrate to the 
T -> limit of Eqs. Q and ©. 

We now show when and how the non-uniform solutions can be found. We look 
for a finite subgroup G of G c , which has to satisfy two simple conditions. We index 
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its irreducible representations by T and use the symbol |T| to specify their dimension. 
Let Pr denote the projectors 

P r = jSE*r(5)5, (H) 
1 1 geG 

with the following properties 

p r p TI = s rTh J2 Pt=l ( 12 ) 

r 

Above, \G\ denotes the cardinal of G and xr(g) the character of g in the representation 
r. Let TLf denote the eigenspace of Hq corresponding to e F . This space is 
invariant to G c and it decomposes according to the irreducible representations of 
G: TLf = ©iPr^^. The sum goes only over those T for which PtTLf ^ 0. In 
general, dimPr^F = «i|r.j|, where ni is the number of representations Tt on TLf- 
The subgroup G we look for must satisfy the following: 

A5) dim Pr^F = |r,|, i.e. we have irreducible representations of G in each Py^TLf- 
A6) \Ti\ = No for some i (we rearrange so that i = 0). 

Let now uq(x) be the particle density when we populate all the states of Hq below 
e F plus the Nq states in Py TLf- Since we assumed that dim7^ > Nq, uq(x) is not 
uniform. The last condition is on the effective potential: 

A7) If Vr ; = (<f>p- 1 Vn 'Pri ) w hh </>p. any norm one vector from Pt/Hf, then 

V Ti - Vr > 0, for alH > 0. (13) 

If the subgroup G satisfying A5-A6 exists and the effective potential satisfies A7, 
then, at least for small A, the zero temperature KS equations have a non-uniform 
solution, 

n(x) =n (x) + o(A). (14) 

This is our main result. 

We end our long introduction with a discussion of the conditions A5-A7. The 
assumption A5 greatly simplifies our proof, but it is not essential (though we don't 
have a proof without A5). A6 assures the closed shell condition and is essential. It can 
be relaxed, for example, we can have two completely filled shells. However, we believe 
that, in the ground state configuration, all No particles occupied the same (lowest) 
energy level. 

Condition A7 refers to the effective potential and it requires, quite naturally, 
that the level populated by the Nq particles to have the lowest energy (in the first 
order in A). We believe A7 is essential. Now, even if we find the correct G, A7 
may not be satisfied, since there is a competition between the Hartree and exchange- 
correlation potentials. For repulsive interactions, the exchange-correlation potential 
must dominate the Hartree potential (see Eq. (1461) ). This is why the electrons 
crystallize at low densities. For attractive interactions (like the Lenard-Jones fermions) 
is viceversa, the Hartree term has to dominate the exchange-correlation. 

To conclude, A5-A6 determines the crystal structure and A7 determines the 
conditions, like the range of densities, in which this structure is stable. 
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Figure 1. If eo and denote the ground state and energy levels below/above e F 
of -ffo, then, at small A, the spectrum of H„ (n 6 S N ) is contained in the hachurcd 
region of the real axis, where: Aeo = A7 a (eo + a), Ae_ = 1 ^2\ 1 ( e - + a )i 
A-e F = \~/a(e° F + a), A+e F = rrg^ (e% + a) and Ae+ = A7 a (e+ + a). We 
define X as the hachurcd region corresponding to some small, fixed A and always 
consider A < A. 



2. The Proof 

The idea behind our proof is the following. We restrict the search for n(x) to the 
densities that are symmetric relative to G, and in a small vicinity of n . Under the 
action of V n , the Fermi level splits into sub- levels, and for A small enough, we show 
that, for all densities in this vicinity, there are exactly N states below e F . This allows 
us to define, for any such n, the density n' corresponding to the potential V ni i.e. a 
map T : n — > n' . The self-consistency means n' — n, i.e. n is a fixed point for T. To 
show that T has a fixed point, we prove that this small vicinity around uq is mapped 
into itself by T and that T is a contraction. 

To define the space L\ ym of G-symmetric densities precisely, we consider the 
isometries 

T g : L\M) -> L\M), (T g n)(x)=n(gx), (15) 

and define 

L\ ym = p| Ker(T g I). (16) 

It is important to notice that L\ ym is a closed subspace of L 1 (A^). Since the solutions 
of the KS equations are not affected if we add a constant to V n , we can assume without 
losing generality that Vr < 0, Vr ; > for i > and min^ = |Vr„|. 

Theorem 2. Let us consider the closed subset of L 1 (A1), 

= {ne Ll ym , \\n\\ L i = N , \\n - n \\ L i < e}. (17) 
Then, for e and A small enough: 

i) The following map is well defined 

: ^i^m — y ^i^m: T e [n] (x) = P^(x,x), (18) 

where P< denotes the spectral projector of H n onto the spectrum below 
(excluding e° F ). 

ii) T c has one and only one fixed point. 

iii) This fixed point is a solution of the KS equations. 
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Proof, i) Let us show first that T e takes S^ y ^ into L l m . Since we exclude e° F , P< 
is well defined for all n £ and can be determined from the resolvent of H n . Also, 

A1-A2 guaranties that the kernel of P< is continuous, thus its diagonal is well defined. 
From A4, H n and consequently P< commutes with all g, g £ G, for all n £ S^™- 
Then 

T e [n]( ff a:) = P<(gx,gx) = (gP<g- V ){x,x) = T e [n](a;). (19) 

Next, we show that T e takes S 1 ^'^ into £ . For this, we need to show that H n 
has exactly N states below e F , for all n £ S^™. For A small, a first, rough location of 
the spectrum can be obtained from Eq. (jSJ. An elementary argument will show that 
the spectrum of H n is always located inside the set I defined and shown in Fig. 1. 
We now investigate the splitting of the Fermi level. For any n £ S^m, the Fermi level 
will split in sub- levels, each corresponding to the different irreducible representations 
Ti (see A5). The energy of any such level can be computed as 

Ev(n) = ^TrP T J z{z-H n )- l ^- v (20) 

with 7' the contour described in Fig. 2. Simple manipulations leads to 

E T {n) =e% + \{<f>°, V n <$) + A 2 /3 r ,A("), (21) 

with (j}^ any norm one vector from PpHp and 

1 F d 

Pr,x(n) = T^TrPr J z(z-i/ )" 1 K(z-if„)- 1 K(z-i/ )- 1 ^.(22) 

7' 

We have an upper bound, Pr,x{ n ) < Pi with /3 independent of T, A or n: 

Pv,\(n) < 2 / NIK* - Ho)-^V n {z - H n )- l V n (z - H^ 1 ^ 

Notice also that 

{<$.,V n <&) = V r + J [V n {x) - V no (x)]\4>° r (x)\ 2 dx. (24) 

Using the eigenvectors expansion of (Hq + a)~ 2 , one can derive 

|^(x)| 2 < (e% + a) 2 k a , (25) 

leading to 

|($,K$)-VH <(e° F + a) 2 k aX e. (26) 

Returning to Eq. it follows from Eqs. and that Er < e° F and E r > e F 
for l^r , as long as 

{e F + ayk a x 

The last thing we need to prove is that, for A small enough, 

IITeW-nolUx <e, (28) 
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Figure 2. The two contours, 7 and 7', surround the states below e^, and the 
states split from the Fermi level, respectively, such that the distance from any 
point on the contours to X is equal to some A > 0. 



for any n in S^i^ and e satisfying Eq. $27). Let us consider 

Z[n} = j( z -H n )-^.+P To J(z-H n )- 1 ^- i , (29) 



denned on the entire S . Based on the previous results, we observe that P< = Z[n] for 

syr, 

if || ||i denotes the trace norm, 



n £ S^ y f n , with e satisfying Eq. l(77|) . Also, notice that n n (x) = Z[n](x,x). Moreover, 



\Z[n] Z[n']\\ x < (1 A °^ )2 II" " n'\\ L i, (30) 



with 



Indeed, let 



7 J l' 



\z + a\ \ 2 \dz\ 



1 + ^ V- ^ 



B = (H + a)- 1 (V n -V n/ )(H Q + a)- 1 , (32) 
and g a ,z{t) = (t + a)/(z — t). After simple manipulations, 

(1 + X(H„ + a)- l V n )- l B(l + \V„,(H„ + a)- 1 )' 1 g^(H n ,) (33) 

and notice that 

\\9aAHn)\\<l+ 1 -^, (34) 

for all n £ S N and z £ 7 or 7'. We can conclude at this step that 

||%]- %']!!!< Aa llflllx. (35) 
(1 - A7J 2 

If we use the polar decomposition AV = S\AV\ of AV = V n — V n >, and define 
A = |AV"| 1 / 2 ( J ff + a) -1 , then from Eq. 

= HAWAII! < ||A*|| ffS ||5A||^ < \\A\\ 2 HS , (36) 

and 

\\A\\ 2 HS = [\AV(x)\(H + a)- 2 (x,x)dx<k aX \\n-n'\\ L1 . (37) 
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With Eq. (|5U|) proven, we can easily end the proof of i). Indeed, for all n £ S^ y f n , 
||T e [n] - n \\ L i = \\Z[n](x,x) - Z[h](x, x)\\ L i 

< \\Z[n] Z[n\h < n Xa f\ 2 \\n-n\\ Ll (38) 

(i - \~f a y 

and 

\\n - 77-H < \\n - n || + ||n - n|| < e + 2JV . (39) 
Thus, Eq. (J2HJ is true if 

Aa X fca(2iVo + e) 

(1 - A 7q )2 - ^ 

and we remark that Eqs. I|27(l and 140(1 can be simultaneously satisfied if A is small 
enough. 

ii) Observe that if we take A small so as to satisfy Eq. (|40ll . then 

Xaxka < 1 (41) 

(1 - A 7 a) 2 ( ' 

Then T e is a contraction, since P< and Z[n] coincide on and 

IITeM-TeMIU^HP^-P^ll! 

= ||Z[n]-Z[n']||i < ||n-nV- ( 42 ) 
Since T e is a contraction on a closed set, it must have one and only one fixed point. 

iii) It follows immediately if we express P£ in terms of the eigenvectors and notice 
that, at a fixed point, n[x) = P^{x, x): 

H n 4> i = c i <t> i , n(x)= \<t>i¥)\ 2 . (43) 

ei<e% 

Together with J n{x)dx = N, the above equations are exactly the KS equations at 
zero temperature. 

3. Examples 

We consider first one of the simplest examples possible: 2N particles on a circle of 
length a. The Kohn-Sham Hamiltonian is H n — —d\ + V n , where x is the coordinate 
along the circle. The ground state of H = —d% is non-degenerate, while all the excited 
states are doubly degenerate. Thus, if we populate the states of H with 27V particles, 
we end up with one particle occupying a double degenerate energy level, containing 
the states a -^ e ikFX and a - 1 / 2 e~ ikFX (e° = k%). 

We now go over the constructions considered in the previous section. The 
continuous group G c are the rotations of the circle and the subgroup G can be taken 
as the identity plus the reflection r : x — > — x. There are two, one dimensional 
irreducible representations of G, x±( r )^ = ^k. The projectors Pr (T — > ±) are simply 
given by 

(P ± f)(x) = -j=[f(x)±f(-x)]. (44) 

They decompose H.f in the invariant, 1 dimensional spaces 

[2 [2 
P+Hf — {\ - cos kpx}, P^TLf = {a/— sin kpx}, (45) 
V a V a 
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Figure 3. a) The molecular structure of the C60 molecule. The small spheres 
represent the carbon atoms. The double bonds are indicated by segments joining 
two atoms, b) A contour plot of the density no(x) defined in the text. The 
dark/lighter regions corresponds to high/low densities. 



each of them providing an irreducible representation for G. Either one of the (±) 
representations can chooses as Tq in A6. We choose the (+) representation, in which 
case, n$(x) = n + a -1 cos(2fci?a;) and condition A7 reads 

J V no (x)cos(2k F x)dx < 0. (46) 

In LDA, if we approximate v xc (no) ~ v xc (n) + v' xc (n)(no — n), Eq. (|46|l reduces to 

v(2k F )+v' xc (n)<0, (47) 

where v is the Fourier transform of the two-body interaction. For QA, Eq. I|46|) simply 
means K(2kp) < 0. We will be led to the same conditions on the effective potential 
if we choose Tq to be the anti-symmetric (-) representation. 

Similar examples can be given for toruscs in higher dimensions. We, however, 
consider the case of fermions on the 2D sphere and show that we can obtain the 
molecular structure of the Cgo molecule. In the Cgo molecule, the carbon atoms sit 
at the points of intersection between an icosahedron and dodecahedron as shown in 
Fig. 3a. There are single and double bonds between the carbon atoms. Since the 
double bond is much stronger than the single bond, we consider C=C as being the 
building blocks of the Ceo molecule. In total, there are 30 double bonds and some of 
them arc shown in Fig. 3a. We consider then 30 point particles (of course, we tried 
60 particles with no success) on a sphere of radius R, described by the Kohn-Sham 
Hamiltonian H n = —R~ 2 L 2 + V n , where Hq = —R~ 2 L 2 is the kinetic energy of a 
particle on the sphere. The energy levels of Hq are simply R~ 2 l(l + 1), / = 0, 1, . . ., 
and if we populate them in order, we end up with 5 particles on the I — 5 level. The 
continuous group is 0(3) and we can take the proper icosahedral group as the finite 
subgroup G. Indeed, under its action, the Fermi level decomposes as 

H F = H u 8 T lu T 2u , (48) 
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i.e. in a 5 dimensional space plus two 3 dimensional spaces (of different symmetries). 
Thus, the proper icosahedral group satisfies A5-A6 with T = H u . Then, if the 
effective potential satisfies A7 (which we expect to happen for certain radiuses R), 
the Kohn-Sham equations for the 30 particles on the sphere have a stable solution 
n(x) = Uq(x) + o(A), where no (a;) is obtained by populating all / < 5 levels plus the 
5 states with I = 5 and H u symmetry. This density is shown in Fig. 3b and the 
resemblance with the Ceo molecule is evident. 

We end this section with a discussion of the numerical results of Ref. 0j on 2D 
electrons and a completely deformable jellium, whose electrostatic potential cancels 
exactly the electrons' Hartree potential. This system is described by the Kohn-Sham 
Hamiltonian: 

H n = -A + v xc (n), (49) 

where v xc (n) is an effective 2D exchange-correlation potential (treated in the local 
density approximation). From the beginning, we should mention that this system 
does not fall in our category. We discussed here systems that are confined in finite 
volumes and the small coupling constant regime. In Eq. (|49(1 . the electrons are free 
to move in the entire 2 dimensional space. The only confining potential is their own 
exchange-correlation potential. Since this potential needs to bind these electrons, we 
cannot talk about the small coupling limit. However, we will show that our analysis 
still applies. We mention that a group analysis of the symmetry breaking in small 
parabolic quantum dots was already carried in Ref. [B]- 

The self-consistent numerical solutions of Eq. (|49|l showed the following: there is a 
first class of stable crystals, which have pure shapes, like triangles, squares and circles, 
and there is a second class with apparently no regularity in the shape. We discuss the 
crystals with pure shape, where the symmetry group is easy to recognize, and the case 
when the electrons are paired and the spin degree of freedom is not important. Clearly, 
for the triangle-shaped crystals, the symmetry has been broken to G — C3V, the group 
that sends an equilateral triangle into itself. For the square ones, G — D^ v , the group 
that sends a square into itself. In Ref. it was found numerically that crystals with 
8 and 10 particles prefer the square geometry. Let us show how one can predict this. 
We first guess a good density n^°\ which is a uniform density inside a square, 
will create a constant, (strongly) negative potential inside the square and the electrons 
will behave, more or less, like they are confined within hard walls. The eigenstates 
for this system are ip n m(x,y) — 4- sin(n7ra;/i) sm(miry/L) (L the size of the square 
and A = L 2 ). The lowest energy state is ipn. The first excited state is double- 
degenerate and corresponds to 1^12 and ip2i- D^ v acts irreducibly on this space. The 
third excited state is non-degenerate and correspond to ^22 • The fourth excited state 
is again double degenerate and corresponds to 1P13 and ip3i : under the representations 
of Di v , this level decomposes into 2 invariant, 1 dimensional spaces (symmetric and 
anti-symmetric combinations of ip\3 and V"3i)- Thus, with 8 particles, we can populate 
the first 3 levels and have the closed shell condition. Thus, the self-consistency can be 
achieved without further reduction of the symmetry. If we write Hq = — A + v xc (n^), 
we reduced the Kohn-Sham Hamiltonian to H n = Hq + V n , where Hq describes, more 
or less, free electrons confined in the square and V n = v xc (n) — v xc (n^) can be 
considered small. The density uq of our theorem is given by populating the first 3 
energy levels of Hq and proving the existence of a self-consistent solution near hq may 
be accomplished now with the methods developed in this paper. The self-consistent 
field will split the fourth energy level, so we can accommodate another pair of electrons 
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without breaking the symmetry, i.e. the crystal with 10 electrons should also have 
square symmetry. Examining either Fig. 2 or the inset of Fig. 3 (where one can see 
the splitting of the second energy level) in Ref. @J, it follows that, in reality, there is 
a slight deviation from the D/± v symmetry. 

4. Discussion 

We can summarize our procedure as follows (hopefully this can be seen in the above 
examples): in the first part, one looks for no, which must be a good approximation 
of the self-consistent density. The main goal is to fulfill the closed shell condition and 
group theoretical methods can be used to accomplish that. After no is found, one can 
construct the map T in a neighborhood of no- In the second part, one investigates if 
T is indeed a contraction near Hq. For small A, the conditions in which this is true are 
already given here. How large this A can be depends on V n and its derivative near no- 

One important remark is that some of these non-uniform solutions can be found 
only if we start the iteration close enough to them. The reason is that the basin of 
attraction of the map T may be small (notice that the usual iterative process of solving 
KS equations consists exactly in constructing the sequence no, T[no], • ■ ■, T m [no], etc.). 
For this reason, there may be additional isomers to the ones found in Ref. 0]. 

As opposed to the Jahn- Teller effect [5] or Pierls instability JUj, which involves 
discrete symmetries and the electronic degeneracies are lifted by a displacement of 
the atoms, the symmetry breaking discussed here involves a continuous group and is 
due solely to the electron-electron interaction, which lifts the electronic degeneracies 
without any change in the external potential. The Jahn- Teller and Pierls instabilities 
may be triggered (but not necessarily) by the instability we discussed here. 

Our analysis does not rule out the existence of more than one self-consistent 
solution of the KS equations. The lowest energy configuration is, of course, associated 
with the ground state, and the higher energy configurations should be associated with 
excited states. In Ref. |7j, the authors find that indeed, since the energy functional 
-E u [n] is not convex, 2?u[n] may have additional extrema, which are excited-state 
densities. 

Degeneracy and symmetry breaking in DFT are well studied concepts in Density 
Functional Theory |111 I12| . and there are many numerical studies on the subject. In 
fact, the modern theory of freezing 1 1 3j . which can be traced back to the pioneering 
work of Ramakrishnan and Yussouff ^1] , is based on the assumption that the liquid- 
solid transitions occurs because of a bifurcation of the same type as we discussed 
here. More exactly, the density of the solid n s is computed self-consistently as the 
linear response of the uniform liquid density m to the introduction of a density change 
An = n s — m. In the current state, the procedure cannot predict the lattice symmetry, 
but rather assumes that the linear response equation has a non-uniform solution, 
with a prescribed crystalline symmetry (usually known from experiment). Applied 
to quantum liquids, this procedure is equivalent to solving the Kohn-Sham equations 
in the quadratic approximation. One can find in Ref. |15| an impressive numerical 
demonstration of the Wigner crystallization of the electron liquid. For finite systems, 
we already mentioned Refs. ^ and [Hj- We add Refs. and |T7) . 

We also want to mention that Bach et al ^Hj have shown that, for any repulsive 
interaction, the energy levels are always fully occupied in the unrestricted Hartree- 
Fock approximation. This result automatically implies that there must be a symmetry 
breaking whenever the last occupied level of H is only partially populated. In 
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contradistinction, Ref. [19] showed that, within the Hartree approximation (ID), there 
is always a symmetry breaking for short, attractive interactions. This definitely shows 
that we have to go beyond the two approximations. 

At the end, we want to mention that we have partial but interesting results in the 
thermodynamic limit and we are currently considering the finite temperature regime. 
We are also looking into the spin case, when the Kohn-Sham equations become self- 
consistent equations for the density n(x) and magnetization vector m(x) . 
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